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Abstract. The global behavior of solutions is described for the 
smallest chemical reaction system that exhibits a Hopf bifurcation, 
discovered in j^. This three-dimensional system is a competitive 
system and a monotone cyclic feedback system. The Poincare- 
Bendixson theory extends to such systems [6l [H |3l [10] and a 
Bendixson criterion exists to rule out periodic orbits [3]. 
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In [7] , Wilhelm and Heinrich discover the "smallest chemical reaction 
system" which may exhibit a Hopf bifurcation. Like the famous Lorenz 
system, it is a three dimensional system with only a single quadratic 
nonlinearity. In [S] , they study the Hopf bifurcation in detail using cen- 
ter manifold and normal form techniques which establish the existence 
of a stable periodic orbit very near the bifurcation point. However, as 
they point out, their analysis does not preclude that chaotic behavior 
may not occur in this simple system. Sprott [12] has cataloged many 
other three dimensional systems with quadratic nonlinearities which 
can have chaotic dynamics. However, we will show that this particu- 
lar chemical reaction system does not support chaos. We give a fairly 
complete analysis of the global dynamics of the system of ordinary dif- 
ferential equations studied in [S] . It is given by the following equations: 



1. Introduction and Main Results 




x' = kx — k2xy 
y' = k^z - ksy 
z' = k^x — k^z 
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where fcj > 0, i = 2, 3, 4, 5 and k = kiA — k4 where A is the concen- 
tration of "the outer reactant of the autocatalytic reaction". The fcj 
are rate constants and k need not be positive; see [7]. Variables x,y,z 
denote concentrations and therefore are nonnegative. 

The scahng x = x/a, y = y/b, z = z/c and the choice k2b = 1, k^c = 
k^b, k^a = k^c leads to the system: 

x' = kx — xy 
(1.2) y' = k^iz-y) 

z' = k^{x — z) 

where we have dropped the bar over variables for simplicity. Our focus 
is on the dynamics exhibited by (11. 2p in the nonnegative octant Mi^. 

System (11. 2p is a competitive system because the signed incidence 
matrix of its Jacobian has a single loop with two positive feedbacks 
and one negative one: 

X ^ z ^ y ^ X 

See [iniE]; the change of variables z —z gives the system the canoni- 
cal form with all off-diagonal entries of the Jacobian being non-positive. 
Since it is three dimensional, the celebrated Poincare-Bendixon The- 
orem extends to its solutions by a result of M.W. Hirsch |2l [TOl [3]: 
A compact omega or alpha limit set that contains no equilibrium is 
a periodic orbit. In fact, system (II. 2p is also a monotone cyclic feed- 
back system |6|, for which the Poincare-Bendixon Theorem also holds 
(regardless of the dimension of the system). 

We note that is the unique equilibrium of (II. 2p if /c < 0; it is 
asymptotically stable if < and unstable when > 0. In the latter 
case, the stable manifold of in Mi^ is 5* = {{x,y,x) G Mi|_ : a; = 0}. 
If > 0, there is an additional equilibrium E = k{l, 1, 1) which is 
asymptotically stable if k < k^ + k^ and unstable if k > k^ + k^ with a 
pair of complex conjugate eigenvalues with positive real part and one 
negative eigenvalue. 

Part (a) of the following result was proved in [8] using a linear Lya- 
punov function. 

Theorem 1.1. The following hold for (II. 2p .- 

(a) If k <0, every solution converges to 0. 

(b) // < k < ks + k^, then every solution except those on S 
converge to E. 

(c) If k > k3 + k^, then every solution not starting on S or on the 
one- dimensional stable manifold of E, converges to a nontrivial 
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periodic orbit. At least one periodic orbit is orhitally asymptot- 
ically stable and the number of periodic orbits is finite. 

Moreover, if k > then any solution not starting in S satisfies 

1 r 

(1.3) lim — / u{s)ds = k, u = x,y,z. 

The main open problem concerning the dynamics of (11. 2p is the 
number of periodic orbits in case (c). Numerical simulations suggest 
there is only one. 

The geometry of the stable manifold of E is described in the following 
result. It is tangent at E to the eigenvector with sign pattern (+, +, — ) 
and has certain monotonicity properties due to the monotonicity of the 
time-reversed system. 

Proposition 1.2. Letk > k^+k^. Then the stable manifold ofE inM.^ 
consists of E and two monotone solutions puif) = {xuit),yuit), Zuit)), t G 
[0,cx)) andpi{t) = {xi{t),yi{t),zi{t)), te[0,cx)). 

Pu satisfies PuifS) = (^^^(O), 0, Zu(O)) with < Xui^i) < k < -2m(0), 
Pu{oo) = E and Xu,yu o-f^ strictly increasing while z^ is strictly de- 
creasing. Therefore, the graph of Pu belongs to [0, A;] x [0, A;] x [A;, oo). 

pi satisfies pi{0) = {xi{0),yi{0),0) with xi{0),yi{0) > k, pi{oo) = E 
and xi,yi are strictly decreasing while zi is strictly increasing. There- 
fore, the graph of pi belongs to [A;, oo) x [A;, oo) x [0,A;]. 

As a consequence of Proposition II. 2[ the unstable manifold of con- 
nects to a nontrivial periodic solution in case (c) of Theorem 11.11 it 
cannot be a heteroclinic orbit connecting to E. 

Theorem 1.3. System (II. 2p has a non-empty compact invariant set A 
that attracts all bounded subsets o/M^. 
Ifk<0, then A = {0}. 

Ifk>0. Then A C [0, Mf where M = MM|M±M. 

Theorem 1.4. Ifk > 0, then ( \1.2\\ is uniformly persistent, i.e., 3t] > 
such that: 

(1.4) min{x(0), ?/(0), z{0)} > ^ liminf mm{x {t),y{t), z{t)} > r]. 

t—^oo 

The attractor A is the disjoint union A = {0} U C U Ai where: 

(i) {0} attracts all bounded subsets of S. 

(ii) Ai is a compact invariant subset of the interior ofM.^ that at- 
tracts all compact subsets ofM.^ that do not intersect S. 

(iii) C is the one- dimensional unstable manifold of {0} which con- 
nects to Ai. 
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2. Proofs 

Proof. Proof of Theorem 11.31 Let (f){t,p) denote the solution of (11. 2p 
which at t = is at p = (x(0), y{0), z{0)). If A; < 0, then x' < and it 
is trivial to show that A = {0}. 

Assume that A; > 0. According to Theorem 2.33 in [11], we must 
show that (11. 2p is point dissipative and for every bounded set B there 
exists tr >0 such that (p{[tr, oo) x B) is bounded. 

We start by showing that it is point dissipative. If x(0) = 0, then 
x{t) = for all t and there exists ri,a > such that y{t) + z{t) < 
ri{y{0) + 2;(0))e~"*. Assume that x(0) > 0. Then direct computation 
shows that 

(2.1) (y/x)' > h{z/x)-{k + h){y/x) 
(z/x)' > h-{k + h){z/x) 

It follows that {z/x)oc > ^ and {y/x)^ > where = 

liminfj^oo fit)- Consequently, for every e > 0, there exists to > such 
that 

X = kx — x^{y/x) < kx — meX^, t > to 

where = (fc^fc^)^^^^^) - e. Hence, x°° = limsup^^^ a;(t) < k/rrie 
and, since e > is arbitrary, x'^ < M. Using this estimate and the 
differential equation for z immediately yields z'^ < M; using this last 
estimate and the differential equation for y yields y°° < M. Thus, we 
have shown that all points are attracted to the bounded set [0, M]^ and 
(11.21) is point dissipative. 

Now we identify a family of positively invariant bounded sets. Let 
< a < < p < and > k/p. Define 

B(a,p,K) = {{x,y,x) e [0,Kf : X = or z/x > a and y/x> p} 

We claim that Bi^c^^p^K) is positively invariant. Indeed, integrating the 
second inequality (12. ip leads to 

(2.2) {z/x){t) > 77^(1 - e-(''+^^)*) + e-('=+'^^)*(z/x)(0) > a 

k ~\~ kc^ 

for any solution starting in B(^„ p K) with x(0) ^ 0. Similarly, using that 
(z/x)(t) > cr in the first inequality (12. ip leads to 

{y/x){t) > ^^(1 - e-^'+'-^^') + e-^'+'^'>\y/x){0) > p. 

As x' = kx — x'^{y/x) < x{k — px) < when x = K, z' = k^{x — 2;) < 
when z = K and < x < K, and y' = k^{z — y) <Q when y = K and 
< z < K, the positive invariance of -B(o-,p,_ft:) follows. 
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Finally, we show that VL > there exists < cr < < p < 

^2^, and K > k/p such that 0(2, [0,Lf) C B(^„^p^K)- Since is 
bounded and positively invariant, this proves that 0([2, oo) x [0, Vf) is 
bounded. Hence, bounded sets have uniformly bounded orbits. 

Let (x(0), 1/(0), 2;(0)) G [0,L]^ and x(0) > 0. Integrating the second 
of (EH) gives 

z{t) ^ k 



>^^(l-e-(^+^^)), t>l. 
Let 0" be the right side of the above inequality and note that a < 



k+ks • 

Substituting this estimate into the first of (12. ip and integrating leads 
to 

x{2)-k + k3^ ' 

Let p be the right hand side of the above inequality and note that 
p < Now we can choose K > k/p so large that x{2),y(2), z{2) < 

K holds for every point (a;(0), y(0), z(0)) G [0,L]^. It follows that 
for every point (a;(0), ?/(0), z(0)) G [0,L]3, we have (x(2), ?/(2), z(2)) G 

B{a,p,K)- n 

Proof. Proof of Proposition II. 2[ The time reversed system of (12. ip is a 
monotone dynamical system with respect to the partial order induced 
by the octant K = {{x,y,z) : x,y > 0, z < 0}. So we make use of 
Theorem 2.8 in [9]. By irreducibility of the Jacobian matrix at E, 
there is an eigenvector v ^ K with nonzero components. Moreover, 
the unstable manifold of E for the time-reversed system consists of E 
and two strictly monotone trajectories Pu & E ~ K and Pi & E + K 
that are described in Theorem 2.8 in [9]. The monotone trajectory P„ 
must exit M5|_ through the y = hyperplane because the derivative of 
y is negative and bounded away from zero. □ 

Proof. Proof of Theorem 11.41 We employ Theorem 8.17 in [11], using 
the notation developed there. Let p = (x, y, z) G Mi]_ and p{p) = 
m.m{x,y, z}. Xq = {p E Mij_ : p{(j)(t,p)) = 0, t > 0} consists of the 
stable manifold 5* of 0. The requisite compactness assumption (H) of 
Theorem 8.17 follows from Theorem 11.31 Q = Upi^xo^ip) = {0} where 
uj{p) denotes the omega limit set of p. M = {0} is obviously an acyclic 
covering of Q and it is an isolated invariant set in Ri]_ because is a 
hyperbolic equilibrium (Hartman-Grobman Theorem). Finally, M is 
weakly p-repelling, meaning that there is no p G M'^ such that p{p) > 
and (f){t,p) — )■ M as t — )■ oo. This follows because the stable manifold 
of is exactly the x = facet S where p = 0. Consequently, we may 
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apply Theorem 8.17 in [H] to conclude that is uniformly weakly p- 
persistent: 3// > such that p{p) > ^ limsup^^^o Pi4>it^p)) > V- We 
conclude uniform strong persistence, i.e., replacing limsup by liminf in 
the above implication, by invoking Theorem 4.5 in [TT]. 

The partition of the attractor A = {0} U C U Ai is a consequence of 
Theorem 5.7 in [11]. □ 

Proof. Proof of Theorem 11.11 Part (a) was proved in [S]. Part (c) 
follows from Theorems 1.1 and 1.2 of p!3] where we check that the 
determinant of the Jacobian at E is —kk^k^ < 0. 

The final assertion uses the uniform persistence established in The- 
orem [LH which implies that the set D = [r/, M]^ is an absorbing set 
for positive solutions of (11. 2p . Since Jq ~ir(^^ = |t ln(x(T)/x(0)) — )■ 
as T — )■ cxD because x is bounded and bounded away from zero, we 
see from the differential equation that MmT^ao ^ loi^ — y{s))ds = 0. 
Indeed, this limit is uniform for all solutions starting in D because 
I ln(a;(T)/a;(0))| is uniformly bounded on D. Thus, 



uniformly for solutions in D . The other limits follow more easily and 
they too are uniform in D. 

Part (b) is more technical. We use the geometric approach of ruling 
out periodic orbits developed in [1] which uses the second additive 
compound Df{p)^'^^ of the Jacobian matrix Df(p), where f{p) denotes 
the vector field fOD : 



See, e.g., formula for Df{p)^'^^ in appendix of [5j. k < k^^ + k^, we 
may find e G (0, 1/2) such that k/{l — 2e) < /cs + k^. Consider the 3x3 
matrix-valued function p — > A{p) where A{p) is a diagonal matrix with 
diagonal entries (1 — 2e)/x, (1 — e)/x, —l/k^. Using the notation of 
[1], we compute the matrix function B(p) = AfA~^ + ADf{py^^A~^ 
where ^/ is the directional derivative of A with respect to vector field 
/. We find that 





h 

k-y -h 
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The Lozinskii measure /i(-B) relative to the norm \p\ = max{|x|, \y\, \z\} 
is (see pg 41 of P): 

X 

fi{B) = max{-eA;3/(l - e), -eh, ^ _ " (^3 + h)} 

Now we use that hmT->oo Jo ^{^)(^^ = k is uniform for all orbits 
starting in the absorbing set B and that k/{l — 2e) < + to 
conclude that 

1 r 

g2 = lim sup sup — / /i(i?((/)(s,p)))c?s < 0. 

t^co peD -L Jo 

Theorem 3.1 of [1] implies that there can be no nontrivial periodic orbit 
of (11.21) . By the Poincare-Bendixson Theorem for competitive three 
dimensional systems [21 [Sj [10] , E attracts all solutions not starting on 
S. □ 
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